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We present 1.1 mm observations for a sample of 16 powerful radio galaxies at 
0.5<z<5.2 and a radio quiet quasar at z=6.3, obtained using the AzTEC bolometer 
array mounted on the ASTE or the JCMT. This paper more than doubles the number 
of high-z radio galaxies imaged at millimetre/sub-millimetre wavelengths. We detect 
probable millimetre- wave counterparts for 11 of the active galaxies. The 6 active galax- 
ies which do not have a probable millimetre counterpart in our images nevertheless 
£f) • have one or more likely associated millimetric source. Thus, we conclude that powerful 

(radio-loud) active galaxies at high-z are beacons for finding luminous millimetre/sub- 
millimetre galaxies at high-z. The flux densities of our AzTEC counterparts imply star 
formation rates ranging from <200 to ^1300 M Q yr _1 . In addition, we find that for 
the radio galaxoes the 1.1 mm flux density is anticorrelated with the largest angular 
size of the radio source. 

We also present new Spitzer imaging observations of several active galaxies in our 
sample. Combining these with archival data, we examine the mid-infrared colours of 
our sample. We find that radio galaxies for which we have detected a probable 1.1 
mm counterpart have mid-infrared colours consistent with dusty starbursts, and are 
usually bluer than high-z Spitzer-sclcctcd active galaxies. In addition, we find arcs of 
24 /im sources extending across ~200-500 kpc, apparently associated with three of the 
radio galaxies. 

Key words: galaxies: high-redshift - galaxies: active - quasars: general - galaxies: 
starburst, sub-millimetre - cosmology: observations - submillimetre - dust 



1 INTRODUCTION vide unique opportunities to study their host galaxy and 

_ . , , . , ,. , . . ^ . _ TT _,_ . . environment (e.g. Venemans et al. 2007). 
Powerful, high-z radio galaxies (z^O.o: HzRGs) continue to 

play a key role in cosmological investigations. Their high , , TT -r,^ , ,, , ., , . , , „ . . , 

^ j j a a a Many HzRGs are embedded withm haloes of ionized 

luminosities at radio and optical wavelengths make them , . , , , ,. , , . . rl 

° gas which are extended on spatiaf scales of tens of kpc, 

useful as beacons for finding massive elliptical galaxies and , , . , ., , . . ,. c 

° i anc i which emit luminous emission lines from various species 

their progenitors out to high redshifts (e.g. Dunlop et al. (<j McCarthy] Spinrad & van Breugel 1995; Reuland et aL 

1996; Rottgering et al. 1997; McLure et al. 1999), and pro- 20Q3a; ViUar _ Martfn et al 2003) Spatially resolved kine _ 

matic studies of these haloes suggest that the cold gas com- 
prising many of these haloes is in infall towards the centre 

* E-mail: ajh@inaoep.mx of the host galaxy (Humphrey et al. 2007), while in others 
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Figure 1. Postage stamps of 1.1 mm AzTEC signal-to-noise maps for the 16 radio galaxies, which show a 160" X 160"field (or 0.9 Mpc 
X 0.9 Mpc to 1.4 Mpc X 1.4 Mpc in physical units) centred on the radio-optical position (+). Both axes are labeled in arcsec, with the 
spatial zero corresponding to the radio-optical position of the radio galaxy. The galaxies are ordered with z increasing left to right, top 
to bottom. Contours start at 2a and increase linearly by la. The FWHM of the beam is shown in the lower left corner of each map. 
Note that the map of 4C+41.17 has a smaller beam FWHM than the rest of the sample, because the observations were made using the 
15 m JCMT (beam FWHM=18"), rather than the 10 m ASTE (beam FWHM=28"). Probable 1.1 mm counterparts for 10 of the radio 
galaxies are detected at ^3<r significance. 



it is being strongly disturbed by the radio jets (e.g. Best, 
Rottgering & Longair 2000; Nesvadba et al. 2006). In most 
cases, ratios between emission lines imply that the radiation 
field of the AGN is the dominant ionization mechanism for 
the haloes, although shocks may make a fractional contri- 
bution in some cases (e.g. Vernet et al. 2001; Humphrey et 



al. 2008). Modelling of the flux ratio Lya/HeIIA1640 mea- 
sured from one-dimensional spectra has revealed an excess 
of Lycv emission above the predictions of AGN photoioniza- 
tion models (Villar-Martm et al. 2007), which suggests the 
presence of star forming regions in the giant gaseous haloes 
(see also Hatch et al. 2009). 



3DSSJ103D+0524 




Figure 2. Postage stamp of the 1.1 mm AzTEC signal-to-noisc 
map for the z=6.3 radio quiet quasar SDSS J1030+0524. The 
160" X 160"field (or 0.9 Mpc X 0.9 Mpc in physical units) is cen- 
tred on the optical position, which is marked by a cross. Contours 
start at 2cr and increase linearly by la. The FWHM of the beam 
is shown in the lower left corner. We do not detect the 1.1 mm 
counterpart of this galaxy. 



Millimetre emission provides an alternative means by 
which the star formation history of HzRGs can be probed. 
Several tens of HzRGs have now been detected at millime- 
tre wavelengths, with the detection rate rising from ~15 
per cent at z<2.5 to >75 per cent at z>2.5 (Archibald 
et al. 2001; Reuland et al. 2003). Millimetric observations 
are sensitive to cold dust that is re-radiating emission re- 
ceived from young stars: implied rates of star formation 
range from a few hundred to ~1500 Mq yr -1 (Dunlop et 
al. 1994; Archibald et al. 2001; Reuland et al. 2004). These 
millimetre/sub-millimetre detected radio galaxies are not 
typical of the wider population of millimetre/sub-millimetre 
selected galaxies (MMGs hereinafter): unlike the vast ma- 
jority of MMGs, the mm emission from radio galaxies is 
often extended on spatial scales of >100 kpc (e.g. Stevens 
et al. 2003), can show morphological complexity (Greve et 
al. 2007), and is not always spatially coincident with the 
radio-optical nucleus of the AGN and host galaxy (Ivison 
et al. 2008). Several HzRGs also appear to have MMGs as 
companions, or are associated with an overdensity of MMGs 
(Ivison et al. 2000; Stevens et al. 2003; see also Priddey, Ivi- 
son & Isaak 2008 and Stevens et al. 2010). 

In this paper we present images and measure 1.1 mm 
flux densities for a sample of 16 powerful radio galaxies at 
0.5<z<5.2 and a radio-quiet quasar at z=6.3, observed with 
the AzTEC bolometer array instrument (Wilson et al. 2008) 
mounted on the Atacama Submillimetre Telescope Experi- 
ment (ASTE: Ezawa et al. 2004) or the James Clerk Maxwell 
Telescope (JCMT), as part of the AzTEC/ ASTE Cluster 
Environment Survey (ACES: Zeballos et al., in preparation). 
Throughout this paper we adopt a flat Universe with Ho=71 
km s -1 Mpc -1 , £Im=0.27 and Aa=0.73. In this cosmology 
l"corresponds to 5.7-8.6 kpc for 0.5<z<6.3. 



2 DATA 

2.1 Sample 

For this investigation we have selected 16 powerful radio 
galaxies at 0.5<z<5.2 from compendia of high-z radio galax- 
ies (e.g. Miley & De Breuck 2008). This sample is designed to 
cover a large range in z, to afford investigations into possible 
cosmic evolution of the properties of powerful radio galaxies 



and their environments. The lower limit to z was set at 0.5 in 
order to minimise the possible contribution to the observer- 
frame 1.1 mm emission from the Sunyaev-Zel'dovich effect 
(SZE), though it is not impossible that the SZE still makes 
a significant (i.e., a few mjy) contribution to the 1.1 mm 
flux densities of radio galaxies at 0.5<z<l (e.g. High et al. 
2010; Menanteau et al. 2010; Foley et al. 2011). However, 
one of our galaxies in that redshift range (MRC 2008-068: 
z=0.55) has been imaged in X-rays (XMM-Newton observa- 
tion ID 0502510501), but shows no evidence for the extended 
X-ray emitting gas needed to produce a significant SZE sig- 
nal. Since we have no evidence for any SZE signal in the 
0.5<z<l radio galaxies, we are unable to make a realistic 
correction for possible SZE contamination. 

Our sample also covers a substantial range in the (pro- 
jected) size of the radio source, with diameters ranging be- 
tween 0.2 kpc and 390 kpc, in order to study the possible 
relationship between the growth or containment of the ra- 
dio source and the millimetric properties of the galaxy. The 
sample purposefully contains some radio galaxies that are 
known to reside in clusters or proto-clusters (e.g. PKS 1138- 
262: Pentericci et al. 1997), some radio galaxies which do 
not appear to reside in overdense environments (e.g. MRC 
2048-272: Venemans et al. 2007), as well as radio galaxies 
for which cluster/overdensity analyses have yet to be under- 
taken. 

In the process of selecting our sample, no strong prefer- 
ence was given to any other observational properties of the 
galaxies, but several strong biases necessarily exist due to 
processes involved in the initial identification of the radio 
galaxies. Firstly, the galaxies are all near the top of the ra- 
dio luminosity function for radio galaxies, and the majority 
of them were selected from radio catalogues on the basis 
of their ultra-steep radio spectra (S v oc v a where a <-l: 
e.g. Chambers et al. 1996). In addition, these galaxies all 
have luminous UV-optical continuum emission, as required 
for the identification of the optical counterpart. They also 
all emit luminous UV-optical emission lines from spatially 
extended regions of cold gas, photoionized by the active nu- 
cleus: such conditions are necessary for obtaining a spectro- 
scopic redshift. In summary, our sample is a subset of high-z 
radio galaxies which have very powerful, steep-spectrum ra- 
dio sources, together with luminous UV-optical continuum 
and line emission. 

It is important to note that roughly 30 per cent of high- 
z radio galaxy candidates do not show bright optical/UV 
emission lines (e.g. Miley & De Breuck 2008), meaning that 
samples selected in the above way may not be completely 
representative of high-z radio galaxies. Interestingly, it has 
been suggested that these 'no-line' radio galaxies are very 
dusty, based on millimetre/sub-millimetre detections thereof 
(Reuland et al. 2003b; Reuland 2005). 

Our sample also includes the z=6.28 radio quiet quasar 
SDSS J103027.10+052455.0 (Fan et al. 2001). Throughout 
this paper, we abbreviate this source name to J1030+0524. 
When the observations were being prepared, this source held 
the record for the active galaxy with the highest known z. It 
is currently unclear whether this quasar resides in a rich en- 
vironment (Stiavelli et al. 2005; Willott et al. 2005, Priddey, 
Ivison & Isaak 2008; Kim et al. 2009). 
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Table 1. Basic properties of the sample of radio galaxies. Columns are as follow. (1) Name of the active galaxy. (2) Active galaxy 
redshift. (3) Rest-frame wavelength, in units of /im, corresponding to 1.1 mm in the observer frame. (4) and (5) are the RA and Dec of 
the active galaxy. (6) The wavelength regime from which the position of the active galaxy was determined, where R = radio, I = mid-IR 
(3.6 and 4.5 (im), O = optical (R-band) and K = near-IR (K-band). (7) 1.1 mm flux density of the millimetre counterpart to the active 
galaxy. Errors in Si.i arc lcr, and upper limits arc 3cr. (8) The observed offset of the mm counterpart relative to the position of the 
active galaxy. In the absence of a probable counterpart we show in parentheses the offset to the closest mm detection. (9) P-statistic of 
the probable mm counterpart. Again, in the absence of a probable counterpart we show in parentheses the P-statistic of the nearest mm 
detection. (10) Log of the far-IR luminosity of mm counterparts. For non-detections, we give 3cr upper limits. (11) The star formation 
rate implied by the far-IR luminosity. Included at the bottom of the table is the radio-quiet quasar SDSS J1030+0524. 
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MRC 2201-555 


0.51 


730 


22:05:04.83 


-55:17:44.0 


I 


6.1+0.8 


2.5 
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12.86 
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MRC 2008-068 


0.55 
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MRC 2322-052 


1.19 
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<12.60 
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-03:44:46.7 
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2.3+0.6 


2.3 


0.0006 


12.68 


480 


MRC 2048-272 


2.06 


360 
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-27:03:03.7 


I 


2.3+0.8 


11.4 


0.02 


12.68 


180 


MRC 0355-037 


2.15 


350 


03:57:48.06 


-03:34:09.5 


I 


3.5+0.7 


6.6 


0.002 


13.02 


1040 


PKS 1138-262 


2.16 


350 


11:40:48.35 


-26:29:08.6 


R 


<3.6 


(14.3) 


(0.007) 


<12.85 


<710 


4C +23.56 


2.48 


320 
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+23:31:45.1 


R 


<1.5 


(21.4) 


(0.005) 


<12.46 


<290 


MRC 2104-242 


2.49 


320 


21:06:58.27 


-24:05:09.1 
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3.7+0.9 


6.9 


0.002 
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0.005 
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2.2 AzTEC 1.1 mm observations 



In the pilot of this programme, the field of 4C+41.17 was 
mapped utilising unchopped raster scans, using the AzTEC 
bolometer array (Wilson et al. 2008) with the James Clerk 
Maxwell Telescope (JCMT) during December 2005. The rest 
of our sample were observed during May to October 2007 or 
July to December 2008, with AzTEC at the Atacama Sub- 
millimetre Telescope Experiment (ASTE), using a Lissajous 
pattern centred on the targeted active galaxy. The obser- 
vations were made under excellent weather conditions, with 
effective atmospheric opacities for each field falling within 
the range 0.03^r^0.06. Integration times varied between 16 
and 35 h per field. The resulting maps cover areas ranging 
from 170 to 300 arcmin 2 for a 50 per cent coverage cut. 



The data reduction procedure is given in Scott et al. 
(2008) with modifications discussed by Scott et al. (2010) 
and Downes et al. (2011). The basic steps are (i) to clean 
the raw time-stream data of spikes due to cosmic rays and 
instrumental glitches; (ii) to then clean using principal com- 
ponent analysis; (iii) to calibrate and bin the bolometer sig- 
nals to produce a map for each individual observation, with 
a pixel size of 2"x2"for 4C+41.17 or 3"x3"for the rest of 
the sample; and (iv) to co-add maps and apply a Wiener 
filter. The observations, the data reduction and the source 
detection/extraction are described in greater detail by Ze- 
ballos et al. (in preparation). 



2.3 Spitzer Imaging 

The Spitzer (Werner et al. 2004) observations of the ac- 
tive galaxies come from both archival data and observations 
from PID 50610 (PI: M. S. Yun). Fifteen of the galaxies 
were observed with the Multiband Imaging Photometer for 
Spitzer (MIPS: Rieke et al. 2004) at 24 /urn, and fourteen 
were observed with the Infrared Array Camera (IRAC: Fazio 
et al. 2004) at 3.6, 4.5, 5.8, and 8.0 /im. Our new observa- 
tions (PID 50610) consist of 5 MIPS maps (MRC 2201-555; 
MRC 2008-068; MRC 0355-037; TN J2009-3040; TN J1338- 
1942) and 3 IRAC maps (MRC 2201-555; MRC 2322-052; 
MRC 0355-037) of HzRG fields. The new IRAC maps arc 
composed of a 4x3 grid of pointing positions, and at each 
position there are 15 dithered frames of 100 sec exposure 
time each. The new MIPS maps are comprised of a 6x6 grid 
of pointing positions, and with each position mapped twice 
witn an individual 30 sec exposure time. 

From the Spitzer Science Center's basic calibrated data 
(BCD) frames we build mosaics for each field using a cus- 
tomized IDL package (Gutermuth et al. 2008). A num- 
ber of common bright source artifacts (banding, jailbars, 
muxbleed, and pulldown effects) are corrected in the indi- 
vidual BCD frames as a first step. Then the code performs 
a rejection of transients and accounts for distortions intro- 
duced by rotation and sub-pixel offsets in building the mo- 
saics. One additional calibration step, performed only on the 
MIPS observations, is a self-calibration to remove systematic 
artifacts (like bright and dark latents) from the maps. For 
each MIPS astronomical observation request ( AOR) , a flat- 
field map is made by taking the median of the background 



at each BCD pixel. Then the median background map is 
normalized, and all BCD images for that AOR are divided 
by the normalized flat field map. The final calibrated, cor- 
rected mosaics for IRAC and MIPS were then resampled to 
a scale of 0.865"pix~ and 1.80"pia; - , respectively. 

Photometric measurements were performed on the mo- 
saics using the IDL routine APER with an aperture diameter 
of 3.8". We apply aperture corrections in the IRAC bands 
of [1.40, 1.38, 1.55, 1.70], and for MIPS we use a correction 
of 2.00. The corrections are determined by selecting isolated 
sources, having no neighboring detection at 5a within 10" 
identified by Source Extractor (Bertin & Arnouts 1996), and 
taking the ratio of the fluxes of those isolated sources with 
apertures of 12.4 and 3.8" diameter. The differences between 
our aperture corrections and those listed in the IRAC In- 
strument Handbook are as one would expected for our 3.8" 
diameter photometry aperture. For each map the median 
of these flux ratios is taken, and then the overall aperture 
correction is chosen as the mean of the set of corrections 
derived for each band. Small field-to-field variations in the 
median aperture corrections introduce some systematic un- 
certainties in our photometry, and so the standard deviation 
of the aperture corrections derived in our maps is added 
in quadrature to the statistical photometric errors for the 
Spitzer bands. In the event that we do not detect an active 
galaxy, we give a 3a upper limit. 

Ten of the radio galaxies in our sample have Spitzer 
photometry reported in the literature (Villar-Martfn et al. 
2006; Seymour et al. 2007). For all but one of these (PKS 
0529-549), our photometry is consistent with the previously 
published measurements. In the case of PKS 0529-549, we 
measure the 24 /im flux to be significantly lower than that 
reported by Seymour et al. (2007): we obtain 632±27 /iJy 
compared to the 942±71 /^Jy of Seymour et al. (2007). The 
higher flux density measured by Seymour et al. is prob- 
ably due to the presence of a second 24 /im source lo- 
cated ~6"SW of the radio galaxy, which would have con- 
tributed significantly to the light measured in their relatively 
large 13"diameter aperture, but significantly less so in our 
3. 8" aperture. 



3 RESULTS 

Figs. 1 and 2 show 160" x 160" 'postage stamps' cut out 
from the full-size signal-to-noise maps, with contours start- 
ing at 2<t and increasing linearly by la. Bars indicate the 
orientation of the radio source, where applicable. Most of 
the active galaxies in our sample have been observed at a 
variety of wavelengths and spatial resolutions, meaning that 
there are various different ways by which one can define the 
position of the galaxy. For this study, we will adopt the co- 
ordinates that are most likely to mark the position of the 
galactic nucleus. Our preferred position is that of the radio 
core, since it corresponds to the active nucleus itself. For ac- 
tive galaxies for which the radio core has not been detected, 
or which are radio quiet, our second choice is to use the cen- 
troid of the SPITZER IRAC emission, averaged across the 
3.6 and 4.5 jj,m bands. For type 2 (i.e., narrow line) objects 
this emission is expected to trace the evolved stellar popu- 
lation, which is likely to be more dynamically relaxed than 
young stellar populations which, if present, could dominate 



at shorter optical-IR and longer (i.e. MIPS) wavelengths. 
For type 1 (broad line) objects, this emission is likely to be 
dominated by the active nucleus. Finally, in the absence of 
both a radio core detection and IRAC data, we adopt the po- 
sition of the optical or near-IR emission. The postage stamp 
images are centred on our fiducial 'radio-optical' position, 
the coordinates of which are listed in Table 1. 



3.1 Association between millimetre sources and 
radio galaxies 

In this paper we adopt a source detection threshold such that 
the peak S/N is >3.0 and also that the spatial FWHM is 
consistent with, or larger than, the full width at half max- 
imum (FWHM) of the filtered point spread function (i.e., 
FWHMJS34"). We have attempted to identify sources asso- 
ciated with our active galaxies as follows. First, we have cal- 
culated the Poisson probability that a detected mm source 
is merely a chance association, rather than being associated 
with the active galaxy. This is given by 

P = 1 - e~ wr2N (1) 

where r is the angular distance of the mm source from the 
optical/IR/radio position of the active galaxy, and N is the 
surface density of mm sources in blank fields that have 1.1 
mm flux densities greater than or equal to that of the de- 
tected mm source (Downes et al. 1986). For N we adopt the 
number counts from the AzTEC/SHADES survey (Auster- 
mann et al. 2010). We consider a mm source to be associated 
with the active galaxy when the P-statistic is less then 0.05, 
i.e., the null hypothesis of chance association is significant 
at the <5 per cent level. All 17 of the active galaxies in our 
sample have at least one associated 1.1 mm source. 

However, we must apply an additional criterion in order 
to determine whether any of the associated mm sources are 
consistent with being at the position of the active galaxy. 
In other words, are the optical/IR/radio and mm positions 
consistent, taking into account the observational errors in 
the position of the mm source? The la positional error of 
an observed mm source is 



o.6 fwhm Y 

\\ y s/N d J 



+ CT 



■p-n I: 



(2) 



where FWHM is that of the filtered beam, S/Nd is the de- 
boosted signal to noise ratio of the millimetre counterpart, 
and o-pnt is the la uncertainty in the pointing of the map 
(adapted from Ivison et al. 2007). We conservatively assume 
that a pnt ~l"(Scott et al. 2010). As an illustration, a source 
with S/Nd=3.0 would have a positional error a pos ~6.9". 
Note that equation (2) does not take into account any sys- 
tematic absolute offset in the map, i.e., residuals after point- 
ing corrections. 

We classify an AzTEC source as the 1.1 mm counter- 
part to the radio galaxy or quasar when P<0.05 and when 
the spatial offset from the radio-optical position r is less 
than 2.7a po £j. There are 11 AzTEC sources which satisfy 
our criteria for classification as 1.1 mm counterparts to the 



1 The positional error distribution is non-Gaussian. 1.1a vos cor- 
responds to a 95 per cent cumulative probability. 
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radio galaxies or quasar. Notice that 6 active galaxies have 
both P^0.05 and r>2.7o pos : this means they have a statisti- 
cally significant association with a millimetre galaxy which 
is not likely to be the active galaxy host. The results from 
this analysis are shown in Table 1. 

Finally, while the presence of other (bright) sources 
within ~30-60"can potentially skew the centroid of a mil- 
limetre source due to the negative side- lobes of the telescope 
beam, our tests show this effect is negligible for our entire 
sample: all such centroid shifts are $Cl"and do not affect our 
counterpart vs. non-counterpart classification. 

3.2 Consistency with previous 
millimetre/sub-millimetre measurements 

A number of the galaxies in our sample have been im- 
aged previously at millimetre or sub-millimetre wavelengths. 
These are 4C+4L17 (Ivison et al. 2000; Stevens et al. 2003; 
Greve et al. 2007), PKS ff 38-262 (Stevens et al. 2003), TN 
J1338-1942 (De Breuck et al. 2004) and SDSS J1030+0524 
(Priddey, Ivison & fsaak 2007). Taking into account the 
(usually) lower spatial resolution of our new observations, 
we find good agreement with the images and photometry 
from the literature. In two cases (PKS 1138-262 and TN 
J1338-1942), several discrete sources detected by SCUBA or 
MAMBO are blended together in our images, due to their 
relatively lower spatial resolution. 

In addition, eight of our sample have been targeted 
with SCUBA photometry (4C+23.56, 4C+41.17: Archibald 
et al. 2001; PKS 1138-262, MRC 0316-257, TN J2007-1316, 
TN J1338-1942, TN J0924-2201: Reuland et al. 2004; SDSS 
1030+0524: Priddey et al. 2003). For seven of these galaxies, 
we find good agreement between our flux density measure- 
ments and the SCUBA photometry. (We use a scaling factor 
of 2.5 between 850^tm and 1.1 mm, assuming Tdust—'iO K 
and /3=1.5: see §3.5.) 

TN J0924-2201 is the exception. Reuland et al. (2004) 
derived a 3<r limit of ^3.2 mjy at 850/^m, which implies a 
1.1 mm flux density of <1.3 mjy, whereas we measure a 1.1 
mm flux density of 3.6±0.9 mjy. The reason for this discrep- 
ancy is not clear. It may suggest that the 1.1 mm source we 
identify with the radio galaxy is not the true 'counterpart', 
and instead is a companion MMG offset from the position 
of the radio galaxy host (see §6). Maps with higher spatial 
resolution and S/N will be needed to address this issue. 

3.3 Stacking 

Next we derive average values of Si.i for various subsets of 
our sample. To do this we stack the AzTEC images at the 
position of the active galaxy, allowing us to include infor- 
mation from the individual observations wherein the active 
galaxy counterpart was not detected at ^3<r significance. 
First, we spatially registered the images. In the case of 4C 
+41.17, we resampled the image to a scale of 3"pixel _1 . 
Next, we assigned to each image a weighting proportional 
to a~ 2 . We then stacked all the images in our sample (Fig. 
3, top left). From this stacked image, we measure an average 
Si. i=3. 4+0. 2 mjy at the position of the active galaxy. We 
repeat this using only the radio galaxy images, and obtain 
an average Si. i =3.2+0.2 mjy (Fig. 3, top right). Repeat- 
ing for the 6 active galaxies for which we did not detect a 



Entire sample 




Figure 3. Postage stamps of the stacked AzTEC images. The 
images were shifted such that the radio-optical position of the 
active galaxy is at the centre, and then co-averaged with weight- 
ings proportional to <r~ 2 . Top left: all 17 active galaxies in our 
sample. Top right: all 16 radio galaxies. Middle left: all non- 
detections. Middle right: non-detected radio galaxies. Bottom left: 
radio galaxies stacked with their radio source axes running hori- 
zontally. Contours begin at 2cr, and increase linearly by 2<r. 



millimetre counterpart, we measure an average flux density 
Si. i=2. 0+0. 3 mjy (Fig. 3, middle left). Finally, for the 5 
undetected radio galaxies we stack to obtain an average flux 
density Si.i=2.4±0.3 mjy (Fig. 3, middle right). 

We have also investigated whether the millimetre emis- 
sion is more extended along the radio axis of the radio 
galaxies in our sample. This is motivated in part by the 
fact that the ultraviolet and optical emission from power- 
ful radio galaxies is usually aligned with the radio source 
(e.g. McCarthy et al. 1987; Chambers, Miley & van Breugel 
1987) , and the apparent alignment between the radio source 
and the large scale distribution of MMGs found by Stevens 
et al. (2003) in a sample of 7 radio galaxies at 2.2^z^4.3. 
To this end, the AzTEC images were registered, rotated 
such that the radio axis runs horizontally, and subsequently 
stacked (Fig. 3, bottom left). Along the radio axis in the 
stacked image, we measure FWHM=39.3±0.6", compared 
to FWHM=40.8±0.6"measured orthogonally to the radio 
axis. Thus, we conclude that there is no significant trend 
for the millimetre emission to be more elongated along the 
radio axis than it is along the orthogonal axis, at least at 
the spatial resolution of our observations. 



3.4 Synchrotron emission 

Before using the our 1.1 mm flux densities to estimate the 
infrared luminosity and the implied rate of star formation, 
we must first consider the impact of synchrotron emission 
from the radio source at 1.1 mm. Indeed, previous stud- 
ies have shown that in some high-z radio galaxies, the syn- 
chrotron emission has sufficiently high luminosity and/or a 
sufficiently flat SED at high frequencies to be detectable in 
the mm-wave regime (e.g. Archibald et al. 2001; Vieira et 
al. 2010). 

To quantify the possible contribution at 1.1 mm from 
synchrotron emission S ex t, we have extrapolated a power law 
of the form S v (xv arad from the two highest frequency pho- 
tometric data points of the observed radio spectra separated 
in v by a factor of ^1.2. When more than one measurement 
is available at 'similar' frequencies, which we define as when 
the ratio between frequencies is ^1.2, we have conservatively 
used the measurement that would result in the highest value 
of S ext (e.g. MRC 0355-037). In Fig 4. we show the radio- 
IR spectral energy distributions of the radio galaxies in our 
sample, including our extrapolations to 1.1 mm of the syn- 
chrotron emission. 

In the first instance, we have used radio photometric 
measurements that represent flux densities integrated across 
the entire spatial extent of the radio source, as these data 
are readily available for all of the radio galaxies in our sam- 
ple. However, extrapolating the radio lobe SED to the mm 
regime is less than ideal. With the possible exceptions be- 
ing the two smallest radio galaxies (MRC 2008-068 and 
TXS 2322-040), this photometry is dominated by the ex- 
tended radio lobes (see e.g. Carilli et al. 1997; Pentericci et 
al. 2000). Lobe spectra usually steepen towards higher fre- 
quencies (e.g. Muxlow, Pelletier & Roland 1988; Murgia et 
al. 1999); therefore, in most cases this extrapolation is ex- 
pected to over-predict the contribution from the extended 
lobes, that is, S ex t is expected to represent an extreme upper 
limit to the contamination due synchrotron emission from 
the lobes. Table 2 gives S ex t alongside Si.i, a ra d and the ra- 
dio wavelengths used for the extrapolations. The table also 
gives references to the photometric data points plotted in 
Fig. 4. 

All 11 of the radio galaxies which we have detected at 
1.1 mm have Si.i measurements exceeding Sext- This means 
that synchrotron emission from the radio lobes clearly can- 
not be responsible for all of the 1.1 mm emission from these 
radio galaxies. Although for 6 of those 11 radio galaxies S ex t 
is greater than the la uncertainty on Si.i, we refrain from 
drawing the naive conclusion that Si.i in some cases may be 
significantly enhanced by synchrotron emission, because our 
extrapolation of S ex t is based on a very favourable set of as- 
sumptions, and the results thereof represent rather extreme 
upper limits. Indeed, in one particular case (MRC 0316-257) 
S ex t is actually a factor of ~2 higher than the 3a lower limit 
to Si.i. 

Moving into the millimetre regime from the radio, the 
flux density of the lobes dimishes more rapidly than that of 
the core (e.g., Carilli et al. 1997), so that at millimetre wave- 
lengths the core typically dominates the spatially integrated 
synchrotron flux density (see, e.g., Figure 1 of Nesvadba et 
al. 2009; Downcs et al. 1996). Thus, we must also consider 
whether the radio core might contribute significantly to the 



1.1 mm emission. The radio core has been detected in only 4 
of our radio galaxy sample (Chambers et al. 1996; Carilli et 
al. 1997; Pentericci et al. 2000, 2001), and for these galaxies 
we have extrapolated the radio core SED to 1.1 mm using 
a powerlaw. One of these galaxies was detected at 1.1 mm 
(MRC 2104-242), and the value of S ex t for its core is several 
orders of magnitude below our Si.i measurement (and well 
below the Sext extrapolated from the radio lobes). 

In summary, the radio lobes are unable to provide suf- 
ficient synchrotron flux at 1.1 mm to be responsible for the 
flux densities we measure at that wavelength, despite our 
use of a very favourable set of assumptions. Furthermore, it 
seems likely that any enhancement by synchrotron emission 
is smaller than the la uncertainty in our flux density mea- 
surements, i.e., insignificant, in essentially all cases. Thus, 
we prefer not to apply any correction for synchrotron con- 
tamination to our measured values of Si.i, because the ex- 
trapolation of the synchrotron flux density is so uncertain. 
Finally, we comment that the issue of possible synchrotron 
contamination would be resolved with observations at longer 
millimetric wavelengths (e.g. 2-3 mm). 



3.5 Implied star formation rates 

In order to convert our measured flux densities into far-IR lu- 
minosities (Lfir) and star-formation rates (SFR) , we adopt 
a grey-body emission template with an emissivity index /3 
and a dust temperature Tdust- This gives 



Mdust = 



SFR- 



Oobs-L/T 



(1 + z)Kv re8t Bv re 
87rhKis rest / kTdus 



t \J dust ) 
/3+4 



(3) 



-2,,/3 



V h 



L f 



10 10 L© 



M e yr~ 



T(/3 + 4)C(/3 + 4)M dust (4) 



(5) 



where Mdust is the mass of dust, S b s is the observed flux 
density, Dl is the luminosity distance, K Vre3t is the mass 
absorption coefficient of the dust at rest-frame frequency 
Vrest, and B l , rl!St (Tdust) is the Planck function at v re st (e.g. 
Archibald et al. 2001). T is the Gamma function, and £ is 
the Riemann zeta function. A Salpeter initial mass func- 
tion has been assumed. Throughout this paper, we assume 
T dust =40 K, 13=1.5 (after Archibald et al. 2001). We adopt 
«375GHz=0.15 m 2 kg -1 (Hildebrand et al. 1983), and ex- 
trapolate to other frequencies assuming k v oc v® (Chini et 
al. 1986). 

The values we have measured for S\.\ imply star for- 
mation rates ranging from <200 to ~1210 Mq yr _1 . The 
sample averaged Si.i we obtained from our stacking analy- 
sis implies an average SFR~600 Mq yr _1 . 



4 COMPARISON WITH RADIO AND 
OPTICAL PROPERTIES 

4.1 Anticorrelation between Si.i and radio source 
size 

It is interesting to examine what, if any, relationship there 
might be between mm/sub-mm flux density measurements 
and the radio source size. This is motivated by a number 
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Figure 4. The observer-frame radio spectral energy distributions of the radio galaxies in our sample and our extrapolation therefrom 
of extreme upper limits on the 1.1 mm flux density due to synchrotron emission. The red vertical bars show our AzTEC 1.1 mm flux 
density measurement, with the vertical length indicating the range of values permitted by the la uncertainty. Where we have not 
detected the radio galaxy's millimetre counterpart, we instead mark its 3cr upper limit with a red arrow. Black diamonds show flux 
density measurements of the synchrotron emission, integrated over the entire radio source. The solid line is our power-law extrapolation 
of this emission from the high-frequency radio regime to 1.1 mm. Where a galaxy's radio core has been detected at more than one radio 
frequency, its flux density measurements are represented by black squares. Our power law extrapolation of the radio core SED to 1.1 
mm is shown by the dashed black line. References to the radio measurements from the literature are given in Table 2. 



of relatively recent observational results suggesting that the 
presence of intense star-formation in high-z radio-loud galax- 
ies is anticorrelated with the size of the radio source (Best, 
Longair & Rottgering 1996; Willott et al. 2002; Humphrey 
et al. 2006). 

For the purpose of this analysis, we define radio source 
size as the largest projected extent of the radio emission: 
we have measured the largest angular extent from published 
radio images, using high spatial resolution images wherever 
possible. Whenever possible we use high spatial resolution 
images (FWHM<0.5"). For sources which are not well re- 
solved, or which do not have a radio image published in the 
literature, we adopt the largest angular size, or the upper 
limit thereto, listed in the literature. Angular sizes were con- 
verted to physical units of kpc, using the conversion dictated 
by our preferred cosmology. 



For high-z radio galaxies, the published radio images 
tend to be sensitive enough only to detect the relatively 
compact and high surface brightness parts of the radio 
source, i.e., hotspots and the core. Therefore, in classical 
double/triple radio galaxies, the largest angular size typi- 
cally represents the angular distance between the most dis- 
tant hotspot, relative to the nucleus, in each of the two jets. 
It is important to be aware that it is quite possible that 
radio emitting structures too faint to be detected lie yet 
farther from the nucleus (see e.g. Coma A: van Breugel et 
al. 1985; 3C171: Blundell 1996; B3 J2330+3927 at z=3.087: 
Perez- Torres & De Breuck 2005). In this case, our measured 
radio source sizes would underestimate the true size. While 
we clearly cannot ascertain whether the hotspot separations 
are genuinely representative of the total radio source sizes 
in our high-z sample, we find it encouraging that in low- 



Table 2. Extrapolation of a powerlaw from radio to millimetre wavelengths. The columns are organised as follows. (1) Name of the radio 
galaxy; if the radio core has been detected at more than one radio wavelength, it is also included in this table. (2) Redshift z of the radio 
galaxy. (3) Observed flux density of the radio galaxy's 1.1 mm counterpart. (4) The extrapolated 1.1 mm flux density of the integrated 
radio source, or radio core where appropriate. (5) Radio spectral index of the integrated radio source, or the radio core where appropriate. 
(6) and (7) The wavelengths of the two radio photometric data points used in the extrapolation. (8) References to photometric data from 
the literature which we have used in this paper; references have been abbreviated to the first two letters of the first author's surname 
and the last two digits of the publication year; the abbreviations arc appended to their corresponding full references at the end of this 
paper. 
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MRC 2201-555 
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0.8 
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-1.51 


1.4 


1.6 


MRC 2322-052 


1.19 
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0.2 
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6.0 
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0.0007 
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PKS 0529-549 
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6.4+0.7 
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-1.15 


1.6 


3.5 
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4.0 


-0.85 


6.0 


11.1 


TN J2009-3040 
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82.0 
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z powerful radio galaxies the hotspot separation is usually 
representative of the total angular extent (see e.g. Black et 
al. 1992). 

In Fig. 5 we plot Si.i versus the size of the radio emis- 
sion for our sample. We find an L-shaped anticorrelation 
between Si.i and radio size. The three galaxies with the 
highest values of Si.i (^5 mJy: MRC 2201-555; MRC 2008- 
068; PKS 0529-549) have relatively small radio sources (s^22 
kpc), while those with relatively large radio sources (>22 
kpc) all have relatively low Si.i (^5 mJy). To test the sig- 
nificance of this apparent anticorrelation in our AzTEC sam- 
ple, we use Cox's proportional hazard model (Isobe, Feigel- 
son & Nelson 1986), which is able to treat censored data, 
i.e., upper limits. Using the null hypothesis that no corre- 
lation is present between Si.i and radio size, we obtain a 
global x 2 =8.3, which corresponds to a significance level for 
independence of 0.003. Thus, the anticorrelation is highly 
significant. 

An alternative test can be made via Poisson probabil- 
ity theory. In this case we adopt a null hypothesis such that 
the population of radio galaxies with relatively large radio 
sources (>50 kpc) has the same fraction of galaxies with 
Si.1 J?5 mJy as does the population of radio galaxies with 
relatively small radio sources (^50 kpc). In our sample, there 
are 3 radio galaxies with Si.i ^5 mJy and D r ad ^50 kpc; 
there are 3 with Si.i <5 mJy and D ra d ^50 kpc; there are 
10 with iSi.i <5 mJy and D ra d >50 kpc; and there are none 
with Si.i ^5 mJy and D ra d >50 kpc. Under the null hy- 



pothesis, the probability of counting radio galaxies with 
Si.i ^=5 mJy and D ra d >50 kpc is 4.5xl0~ 5 . Therefore, we 
reject the null hypothesis, and conclude that the AzTEC ra- 
dio galaxy sample shows a highly significant anticorrelation 
between the 1.1 mm flux density and the apparent size of 
the radio source. 

Shown in Fig. 6 is Si.i versus radio size for the AzTEC 
counterparts together with HzRGs that have been observed 
at 850^m using SCUBA (Archibald et al. 2001; Reuland et 
al. 2003b, 2004; Stevens et al. 2003). SCUBA 850^m flux 
densities have been scaled by 1/2.5 to extrapolate their 1.1 
mm flux densities. This scaling factor is calculated assum- 
ing T du st='i0 K and 0=1.5 (§3.5). Where an HzRG has been 
observed more than once at 1.1 mm and/or 850/im, we use 
the most sensitive of the flux density measurements. Us- 
ing Cox's proportional hazard model we find this combined 
dataset does not show a significant anticorrelation between 
Si.i or SFR and radio size (P=0.18), in agreement with the 
statistical analysis of Reuland et al. (2004). However, us- 
ing Poisson probability theory we find that the absence of 
HzRGs at Si.i ^2.5 mJy and D ra d >200 kpc is significant 
at the 98 per cent level (P=0.02). 

It seems plausible that this reduction in the significance 
of the trend is due to the different observational methods 
employed in our AzTEC/ASTE survey of HzRGs, compared 
to those predominantly used for the several SCUBA/JCMT 
surveys. Whereas our maps allows us to reject sources that 
are nearby to the radio galaxy's position but unlikely to 
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Figure 5. The anticorrelation between 1.1 mm flux density and 
radio source diameter in our AzTEC sample. The dashed lines 
show the divisions used to calculate the Poisson probability 
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Figure 6. Similar to Fig. 5, but including also HzRGs with mea- 
surements or limits at 850/im. We have extrapolated 850^tm flux 
densities to 1.1 mm by multiplying by 1/2.5. To avoid complicat- 
ing the figure, error bars are not shown. 



be the counterpart, the SCUBA photometry employed by 
Archibald et al. (2001) and Reuland et al. (2004) simply 
obtains the flux density from a region of sky sampled by 
the telescope beam, and emission from a bright nearby com- 
panion MMG may contaminate such measurements, thereby 
adding scatter to the trend. 

Finally, we point out that the trend discussed above 
would not change significantly if we were to include 
only millimetre/sub- millimetre measurements obtained from 
maps (e.g. from Ivison et al. 2000; Stevens et al. 2003; De 
Breuck et al. 2004; Greve et al. 2007; and this paper), i.e., 
excluding measurements obtained from simple SCUBA pho- 
tometry. We would still find that, of the radio galaxies with 
radio sources smaller than 200 kpc, 50 per cent have Si.i <5 
mjy, while the radio galaxies with radio sources larger than 
200 kpc all have Si.i <5 mjy. In §6 we discuss the origin 
and implications of this result. 
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Figure 7. Flux density at 1.1 mm versus Lya/Hell. Dashed ver- 
tical lines separate three different areas on the plot: Lyo/Hell 
lower than predicted by AGN photoionization models (absorbed); 
Lya/Hell consistent with AGN photoionization models; and 
Lya/Hell higher than predicted by AGN photoionization mod- 
els (enhanced). Black diamonds show the high-z radio galaxies 
from our sample which have Lyo/Hell measurements in the lit- 
erature. Red triangles show high-z radio galaxies from the litera- 
ture which have both SCUBA 850^tm flux density measurements, 
which we have scaled by a factor of 0.4 to account for the shorter 
wavelength, and Lyo/Hell measurements, the latter from Villar- 
Martm et al. (2007). 



4.2 The Lya emission nebulae 

In addition to heating the surrounding dust young stars also 
produce HII regions, which emit strong UV-optical emission 
lines. The strongest of these lines is Lya A1216, and its lu- 
minosity is primarily a function of the number of ionizing 
photons absorbed, and the quantity and geometry of dust. 
It is relatively insensitive to the metallicity or the gas den- 
sity of the HII regions. Thus, Lya carries useful information 
about dust associated with star forming regions. 

At z>2, Lya is redshifted into the optical regime, mak- 
ing this line accessible to ground-based telescopes. Indeed, 
all 11 of the z^2.15 radio galaxies in our sample have mea- 
surements in the literature for the narrow (FWHM<3000 
km s _1 ) Lya emission. In table 3, we list the Lya luminos- 
ity for these 11 radio galaxies, measured from narrow band 
images (e.g. Knopp & Chambers 1997; Reuland et al. 2003a; 
Venemans et al. 2007), and/or from spectra with a long-slit 
placed along the major axis of the radio emission (e.g. De 
Breuck et al. 2000a). 

It is interesting to compare these measured Lya lumi- 
nosities against the expected (intrinsic) Lya luminosity of a 
massive starburst. First, we derive an approximate relation- 
ship between the Lya luminosity and SFR or Lfir, as 






W^iSFR/WOM, 

-,44/ r , 1n l2 



@yr 

1O m (L j , / r/1O"L ) ergs 



ergs 

i 



(6) 
(7) 

(8) 



which simplifies to become 
Luyoc/LpiR ~ 0.03 

under the assumptions that (i) the starburst is continuous 
and has an age in the range ~l-5 Myr; (ii) all ionizing pho- 
tons produced by the stars are absorbed in the surround- 
ing HII region; and (iii) 67 per cent of photo-ionization 
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events lead to the emission of a Lya photon (Binette et 
al. 1993). Of the 8 galaxies in our sample which have been 
detected in both Lya emission and in 1.1 mm continuum 
emission, 3 have values of LLya/LFiR which are in order-of- 
magnitude agreement with the expected value ~0.03, and 5 
have LLya/LFiR values more than an order of magnitude 
lower (Table 3). 

Thus far, we have ignored another luminous source of 
ionizing radiation which is likely to contribute significantly 
to the production of Lya photons - the active nucleus. The 
ratio between Lya and Hell A1640 is sensitive to the SED 
of the ionizing continuum and, therefore, it can be used to 
assess the extent to which stellar-photoionized HII regions 
contribute to the UV-optical emission line spectrum. Mod- 
els for pure AGN photoionization predict Lya/Hell = 8 - 
15 (e.g. Humphrey et al. 2008), depending on the spectral 
slope of the ionizing continuum between 13.6 eV and 54.4 
eV. On the other hand, photoionization by young stars does 
not result in significant Hell emission relative to Lya. Thus, 
a Lya/Hell ratio significantly in excess of 15 would suggest 
a significant contribution from stellar-photoionized HII re- 
gions. A number of HzRGs, including 3 of our sample (TN 
J2009-3040, 4C+41.17 and TN J1338-1942), do indeed show 
such an excess of Lya emission. Villar-Martm et al. (2007) 
examined the Lya/Hell and Lya/CIV ratios of 61 radio 
galaxies at 1.79^z^4.41, and identified a Lya excess in 11 
of their sample; they calculated that a SFR of ~200 Mq yr _1 
would be required to explain the most extreme objects. 

Naively, one might expect there to be a correlation be- 
tween detection of a Lya excess and detection of millimetre 
continuum, since both are thought to be produced by pow- 
erful starbursts. In Fig. 7 we plot 5i.i against the Lya/Hell 
ratio for radio galaxies from our sample which have Lya 
and Hell detections. We also plot radio galaxies which have 
been observed at 850/im (from Table 1 of Villar-Martm et 
al. 2007), under the assumption that S'i.i=5'85o/2.5. Con- 
trary to initial expectations, we find no trend between 1S1.1 
and Lya/Hell. The possible origin of this result is discussed 
in 86. 



5 ANALYSIS OF THE MID-IR SPECTRAL 
ENERGY DISTRIBUTION 

In Fig. 8 we show the Spitzer IRAC [3.6]-[4.5] versus [5.8- 
8.0] colour-diagram, adapted from Stern et al. (2005), for 
all HzRGs in our sample that are detected in all 4 IRAC 
channels. Most of our target HzRGs, as well as those studied 
by Seymour et al. (2007), are found within the region of 
red rest-frame optical and near-IR color characterisitic of 
hot dust, overlapping with the distribution of the Spitzer- 
selected power- law QSOs of Lacy et al. (2004) and Martfnez- 
Sansigre et al. (2008). While the colours of our target HzRGs 
are marginally consistent with colours expected for active 
galaxies, those detected by AzTEC (filled squares) appear 
near the edges of the distribution, along the model tracks of 
young dust obscured stellar clusters modeled by Efstathiou 
et al. (2000). This is also the area where the majority of 
submillimetre-bright sources are clustered (Yun et al. 2008, 
2011), and these AzTEC detected HzRGs display rest-frame 
near-IR properties more like SMGs than other HzRGs. 

The diagram shown in Fig. 9 aims to distinguish ac- 
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Figure 8. Spitzer IRAC [3.6]-[4.5] versus [5.8]-[8.0] colour-colour 
diagram for the radio galaxies in our sample which have been 
detected in all 4 bands (adapted from Stern et al. 2005). Solid 
squares arc HzRGs for which we have detected a probable 1.1 
mm counterpart; open squares are HzRGs for which we have not 
detected a probable 1.1 mm counterpart. The small open circles 
are the HzRG sample of Seymour et al. (2007). Asterisks (*) show 
the infrared QSOs with power law spectra from Lacy et al. (2004) 
and Martinez-Sansigre et al. (2008). The dot-dashed line indicates 
the AGN region proposed by Stern et al. (2005). The long-dashed 
black line shows the MMG counterpart region proposed by Yun 
et al. (2008). The dotted line and the solid line show the redshift 
evolution colour-colour track for a starburst age of 37 and 72 Myr, 
respectively, based on the theoretical starburst spectral energy 
distribution models of Efstathiou et al. (2000); crosses (x) along 
the 72 Myr track mark z = 0.5, 1, 2, 3, 4, and 5. The thick 
short-dashed line represents power-law spectrum sources having 
a=0.3-1.0, where S v <x v~ a . 



tive galaxies from starburst systems using two different di- 
agnostic criteria proposed recently by different groups: (1) 
S24/SR > 1000 (Dark Optical Galaxies, or DOGs: Dey et 
al. 2008; Fiore et al. 2008); and (2) Ss.o/S-i.s > 2 (Ivison et 
al. 2004; Pope et al. 2008). The DOGs were first identified 
in the study of a population of z ~ 2 infrared bright galaxies 
that are extremely faint in the optical bands, and Fiore et al. 
(2008) concluded based on their X-ray analysis that as many 
as 80% of these dust-obscured galaxies host a Compton-thick 
AGN. The latter criterion was initially developed to exploit 
the presence of PAH features and a power-law AGN spectral 
shape in the rest-frame near-IR, similar to the information 
used in Fig. 8. 

In this diagram, the majority of the HzRGs detected 
in MIPS 24 [ira by Seymour et al. (2007) appear above 
the DOGs line of S 2 4/S R > 1000 (dotted line), suggest- 
ing warm dust and a heavily obscured AGN are present in 
these HzRGs. All six of our target HzRGs detected in the 
MIPS 24 u.m bands are below the dotted line and are on 
average bluer than the others. In fact, the three AzTEC- 
detected HzRGs are among the bluest HzRGs. The AzTEC 
detected HzRGs also show a flatter S8.0/S4.5 flux ratio than 
the rest of the HzRGs. Taken together, these rest-frame opti- 
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Table 3. Radio and optical properties of the radio galaxies. A dash (-) indicates that a quantity has not yet been measured, or is not 
available in the literature. Columns are as follows. (1) Source name. (2) Source redshift. (3) The 1.1 mm flux density of the millimetre 
counterpart, measured from our AzTEC data. (4) The spectroscopic Lyo luminosity in units of 10 44 erg s _1 . (5) Lya luminosity in units 
of 10 44 erg s — 1 measured from narrow-band images. (6) The flux ratio Lya/Hell (7) Ratio of the Lyo luminosity to the far-infrared 
luminosity; the theoretically expected value is ~0.03; Ly« luminosities measured from narrow band images are used where available, 
otherwise luminosities measured from long-slit spectra are used. (8) Projected size of the radio source (D ra dio)- (9) References for the size 
of the radio source, and for the luminosities of the Lya and Hell emission lines; references have been abbreviated to the first two letters 
of the first author's surname and the last two digits of the publication year, and the abbreviations are appended to the full bibliographic 
references at the end of this paper. 



Source name 


z 


Si.i 


l^Lya 


(slit) 


Lhya (image) 


Lya/HcII 


LLya/LFIR. 


Dradio 


Refs 






(mjy) 


(10 44 


erg s -1 ^ 


) (10 44 erg s" 1 ) 






(kpc) 




(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


MRC 2201-555 


0.51 


6.1+0.8 


- 




- 


- 


- 


22 


Bu06 


MRC 2008-068 


0.55 


8.6+0.9 


- 




- 


- 


- 


0.2 


JeOO 


MRC 2322-052 


1.19 


<2.0 


- 




- 


- 


- 


64 


Be99 


TXS 2322-040 


1.51 


2.3+0.6 


- 




- 


- 


- 


0.5 


Li07 


MRC 2048-272 


2.06 


2.3+0.8 


- 




0.65 


- 


0.0035 


70 


PeOO Ve07 


MRC 0355-037 


2.15 


5.0+0.8 


0.40 




- 


.3.0 


0.0010 


101 


Go05 R697 


PKS 1138-262 


2.16 


<3.4 


1.7 




25 


10.7 


>0.025 


80 


Ca97 R697 Vc07 


4C +23.56 


2.48 


<1.4 


0.40 




2.1 


5.3 


>0.018 


390 


Ca97 Ch96 R697 Ci98 Kn97 


MRC 2104-242 


2.49 


3.7+0.9 


2.9 




- 


13 


0.010 


180 


PeOO Vi99 OvOl 


PKS 0529-549 


2.58 


6.4+0.7 


0.41 




- 


12.3 


0.00084 


10 


Br07 R697 


MRC 0316-257 


3.13 


<1.9 


0.21 




0.7 


- 


>0.0047 


59 


Ca97 At98 Vc07 


TN J2009-3040 


3.16 


3.3+0.9 


1.4 




3.0 


24.3 


0.013 


54 


Bo07 DcBrOO Ve07 


4C +41.17 


3.79 


3.5+1.0 


2.1 




13 


26.5 


0.054 


145 


Ch90 Ch96 R697 Re03 


TN J2007-1316 


3.83 


2.8+0.9 


0.26 




- 


9.0 


0.0014 


26 


Bo07 DeBrOO 


TN J1338-1942 


4.10 


<4.7 


1.7 




4.5 


18.4 


>0.014 


80 


PeOO DeBr99 DeBrOl Vc07 


TN J0924-2201 


5.19 


3.6+0.9 


0.10 




0.15 


- 


0.00074 


8 


vB99 DeBrOl Vc07 



cal and near- and mid-IR diagnostic colours suggest that the 
AzTEC-detected HzRGs have bluer colours of star-forming 
galaxies than those of dust-obscured AGNs with charac- 
teristic hot dust emission commonly associated with other 
HzRGs. We note that the overall rest frame optical colours of 
even the most dusty star forming galaxies, such as ultralumi- 
nous infrared galaxies in the local universe, are surprisingly 
blue because young stellar populations distributed through- 
out the galaxy dominates the overall colour and luminosity 
(see Chen et al. 2010). 

In summary, the HzRGs which have an AzTEC coun- 
terpart also have mid-infrared spectral energy distributions 
that are consistent with dust-obscured starburst systems. 



5.1 24/im images 

In Fig. 10 we show 160" x 160"postage stamps of the 24/im 
emission, with contours of the 1.1 mm emission overlayed, in 
order to illustrate interesting morphological features around 
3 of the radio galaxies. A detailed analysis of clustering in 
the full 170-300 arcmin 2 fields (as opposed to the 7 arcmin 2 
'postage stamps' considered herein) of our sample of 17 high- 
z active galaxies will be presented in a future paper (Zeballos 
et al., in preparation). 



5.1.1 MRC 2201-555 

The 24 lira MIPS image shows an arc of emission running 
through the position of the radio galaxy, with a total spatial 
extent of ~38"or 230 kpc. The NW portion of the arc takes 
the form of a relatively bright source ~ 16" (100 kpc) to the 
NW, which is linked to the radio galaxy by a bridge of lower 



surface brightness emission. This structure is too compact to 
resolve at the resolution of our AzTEC/ASTE observation. 

5.1.2 MRC 0355-037 

The MIPS image shows an arc comprised of several 24 fim 
sources extending ~60"(~500 kpc) to the WSW from the 
position of the radio galaxy. At or near to the end of this 
arc, our AzTEC map shows a peak which has S/N=2.6, 
after correcting for the negative sidelobes of radio galaxy's 
AzTEC counterpart. 

5.1.3 TN J 2001 -1318 

An arc of several 24 fim sources extends ~50"(~360 kpc) 
SSE from the radio galaxy. Similarly, in our AztEC image 
there is arc of 1.1 mm emission, cospatial with the 24 fim 
arc, with peak S/N=2.4 after correcting for the negative 
sidelobes of radio galaxy's 1.1 mm counterpart. 



6 DISCUSSION 

In the process of this study, we have identified a statisti- 
cally significant anti-correlation between the size of the radio 
source and the brightness of the millimetre continuum emis- 
sion (or SFR). This is by no means the first investigation 
into the possible relationships between the size of a radio 
source and the other properties of radio galaxies; numerous 
relationships have been found previously between the ob- 
served size of powerful radio sources and various UV, optical 
and infrared properties. Best, Longair & Rottgering (1996) 
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Table 4. Photometric measurements from the Spitzer observations of our sample of high-z active galaxies. 



Source name 


z 


S3. 6 


S4.5 


S5.8 


S8.0 


S24 






(mJy) 


(/iJy) 


0*Jy) 


(A*Jy) 


0*Jy) 


MRC 2201-555 


0.51 


377+6 


401+7 


522+22 


640+27 


1547+255 


MRC 2322-052 


1.19 


147+2 


124+2 


96+5 


90+5 


- 


TXS 2322-040 


1.51 


- 


- 


- 


- 


51+15 


MRC 2048-272 


2.06 


61+2 


71+3 


91+14 


73+20 


268+61 


MRC 0355-037 


2.15 


32+1 


35+1 


42+4 


48+5 


812+136 


PKS 1138-262 


2.16 


353+5 


522+9 


1048+44 


2024+85 


5111+824 


4C +23.56 


2.48 


60+1 


80+1 


164+7 


505+21 


5409+876 


MRC 2104-242 


2.49 


27+6 


33+7 


37+22 


<75 


876+255 


PKS 0529-549 


2.58 


46+2 


55+2 


68+12 


102+14 


632+110 


MRC 0316-257 


3.13 


27.1+0.5 


27.2+0.6 


30+2 


57+3 


661+118 


TN J2009-3040 


3.16 


- 


- 


- 


- 


18+4 


4C +41.17 


3.79 


21.3+0.6 


26.0+0.7 


36+3 


52+4 


467+82 


TN J2007-1316 


3.83 


62+2 


61+3 


65+14 


172+20 


594+102 


TN J1338-1942 


4.10 


21.8+0.8 


12.8+0.6 


17+3 


17+3 


283+61 


TN J0924-2201 


5.19 


6.7+0.3 


10+3 


9.7+0.9 


<37 


217+68 


SDSS J1030+0524 


6.28 


86+1 


109+2 


92+5 


109+7 


523+92 
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Figure 9. Log S24/SR versus Ss.o/Si.5 for the radio galaxies in 
our sample which have been detected in the Spitzer MIPS 24 (im 
channel, in the IRAC 4.5 and 8.0 fim channels, and in the optical 
R-band. Filled squares show radio galaxies for which we have de- 
tected a probable counterpart at 1.1 mm; open squares show radio 
galaxies for which we have not detected a probably 1.1 mm coun- 
terpart. The small circles are the HzRG sample of Seymour et al. 
(2007) as in Figure 8. Note that Sr is the flux density correspond- 
ing to the R or r magnitude of the radio galaxy. References to the 
optical photometry are as follows. MRC 2201-555: Bu06. MRC 
0355-037: Go05. PKS 1138-262: Pe97 MRC 0316-257: Mc96. 4C 
+41.17: Mi92. TN J1338-1942: DeBr02. Starburst galaxies (SB) 
are expected to lie at Log S24/SR < 3, Sg.o/Si.s <2, whereas 
active galaxies are expected to lie above and/or to the right of 
this region of colour space (see Pope et al., 2008; Dey et al. 2008; 
Fiore et al. 2008). 



found a morphological evolution with radio size in their sam- 
ple of eight 3CR radio galaxies at l<z<1.3. For galaxies 
with relatively smaller radio sources, their HST images re- 
veal a string of bright knots tightly aligned with the radio 
emission. The galaxies with relatively larger radio sources, 
while still showing the alignment effect (e.g. Chambers, Mi- 
ley & van Breugel 1987; McCarthy et al. 1987), are generally 
more compact and contain fewer bright optical knots. Fur- 
thermore, galaxies with smaller radio sources, the extended 
emission line region (EELR) tend to have a relatively smaller 
observed spatial extent, as measured with Lya (van Ojik 
et al. 1997). The EELR of these galaxies also show more 
extreme emission line kinematics, have brighter UV-optical 
emission lines, and ratios between emission lines that imply 
they have lower ionization states (Best, Rottgering & Lon- 
gair 2000; Inskip et al. 2002; Humphrey et al. 2006). In addi- 
tion, the polarization of the rest-frame UV continuum emis- 
sion is correlated with radio source size (Solorzano-Inarrea 
et al. 2004), suggesting that the ratio between scattered light 
and starlight is relatively higher in radio galaxies with larger 
radio sources. The anticorrelation between 850 /im luminos- 
ity and UV continuum polarization found by Reuland et al. 
(2004) appears to confirm that the UV continuum polariza- 
tion is indeed strongly dependent on the rate of star forma- 
tion. Furthermore, HzRGs showing anomalously high values 
of Lya/Hell (>15), likely due to the presence of a powerful 
starburst (~100 Mq yr _1 ) tend to have smaller radio sources 
than those with 'normal' or 'low' values (Villar-Martm et al. 
2007). These radio-optical trends are commonly supposed to 
be the result of (a) interaction between the radio source and 
the host galaxy, (b) slowing of radio source expansion in ra- 
dio galaxies that are in relatively denser environments, or 
(c) age effects of the radio source. 

A previously unanswered question was whether the star 
formation rate is intrinsically higher in HzRGs with smaller 
radio sources, or whether instead the UV-optically derived 
SFR merely appears higher due to some difference in dusti- 
ness, geometry, etc. Observations in the millimetre/sub-mm 
regime give us a new perspective from which we may address 
this question, since the observations are sensitive to UV- 
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Figure 10. Postage stamps of the 24/xm emission (greyscale) 
with contours of the 1.1 mm emission overlayed. For MRC 2201- 
555, the AzTEC contours are the same as in Fig. 1. For MRC 
0355-037, the contours start at la and increase by la. For TN 
J2007-1316, the contours are 1.7, 2.2, 2.7 a. North and East are at 
the top and left of the images, respectively. The cross (+) marks 
the position of the active galaxy. 
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Figure 11. Radio size as a function of z. In the interest of clarity, 
error bars are not shown. Note the clear anticorrclation between 
these two parameters, which complicates interpretation of any 
observed evolution of properties with z (see text). 



optically obscured (rather than unobscured) star forming 
regions. Assuming that the millimetre/sub- mm continuum is 
dominated by thermal emission from starburst-heated dust, 
as we believe to be the case, then our anticorrelation be- 
tween Si.i and radio size leads us to conclude that HzRGs 
with relatively smaller radio sources do indeed tend to have 
intrinsically higher SFR. 

In light of this conclusion, it seems natural to now ques- 
tion whether the increase in the mm/sub-mm detectability 
of HzRGs at z^2.5 is due to a genuine redshift evolution in 
the rate or mode of star formation in the host galaxies, as has 
commonly been supposed (Archibald et al. 2001; Reuland et 
al. 2004). Indeed, we propose an alternative hypothesis: that 
there is a bias towards identifying HzRGs with smaller radio 
sources at higher z (Fig. 12) because these sources are asso- 
ciated with more luminous narrow emission line regions (e.g. 
Best, Rottgering & Longair 2000). And since smaller radio 
sources typically have more luminous mm/sub-mm emission 
(Figs. 5 & 6), we propose that this translates into a bias 
in favour of identifying HzRGs with significantly more lu- 
minous mm/sub-mm emission at higher z. Identifying and 
studying the 'missing', line- weak HzRGs will be crucial for 
testing this hypothesis (see Reuland et al. 2003b; Reuland 
2005; Miley & De Breuck 2008). 

It is not known whether the observed size of the radio 
source is determined primarily by the density of the inter- 
stellar/intergalactic medium through which the source prop- 
agates, or whether the size instead determined primarily by 
the age of the radio source. If we were to assume the lat- 
ter, then we would be able to gauge the lifetime of luminous, 
dusty star formation episodes (Si.i >2.5 mjy; SFR>500 M© 
yr _1 ) closely associated with HzRGs. In Fig. 5, all of the lu- 
minous millimetre counterparts are associated with HzRGs 
that have radio source diameters of d<200 kpc: we adopt 
d-~200 kpc as the approximate stage in the growth of the 
radio source at which the luminous millimetre counterparts 
fade to Si.i <2.5 mjy (SFR<500 M Q yr" 1 ). Assuming a 
conservative expansion speed of ~0.01c for the radio source 
(e.g. Best et al. 1995), this diameter implies an age of ~30 
Myr. A similar age was estimated by Willott et al. (2002) 
for a smaller sample of high-z, radio-loud quasars. However, 
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we must be cautious about generalising this to obscured star 
formation in non-active galaxies (i.e., to sub-mm galaxies in 
general): HzRGs contain two powerful sources of feedback, 
the radio source and the accretion disc, both of which are 
able to quench star formation and, perhaps, to ignite it. 
Concievably, feedback of this kind could either lengthen or 
shorten the lifetimes of star formation activity in HzRGs, 
and could mean that lifetimes of obscured star formation 
activity determined for HzRGs are not necessarily repre- 
sentative of those for millimetre/sub-millimetre galaxies in 
general. 

We have also investigated whether there might be a 
correlation between the millimetre flux density and the 
Lya/Hell ratio, which is expected to be sensitive to SFR 
also. While Lya can originate from stellar photoionized HII 
regions and from AGN-ionized nebulae, the Hell line is ex- 
pected to come predominantly from nebulae that are ionized 
by the AGN. Though both Lya/Hell and the millimetre flux 
density are expected to be correlated with SFR, and there- 
fore with each other, no significant trend was found between 
them. Furthermore, the luminosity of the Lya emission is of- 
ten more than an order of magnitude lower than one would 
have expected given the SFR derived from the 1.1 millimetre 
flux density (we estimate LL ya / 'Lfir ~0. 0007-0. 054, while 
the expected value is ~0.03). This may mean that one (or 
both) of the millimetre and Lya emission is an inaccurate 
tracer of the SFR in HzRGs. This could be due to destruc- 
tion of Lya photons by the dust associated with this kind 
of star formation activity. 

Thus far, we have tacitly assumed that the sources we 
identify as 1.1 mm counterparts are the radio galaxy hosts 
themselves. However, it is important to be aware that the 
spatial resolution of our observations, together with the sig- 
nal to noise ratios of the detections, has resulted in lcr posi- 
tional uncertainties of ~2"-16"for our 1.1 mm counterparts. 
This means that the sources we identify as 1.1 mm counter- 
parts could in actual fact be offset by tens or even hundreds 
of kpc from the galactic nucleus, which would place them 
in the outskirts or outside of the host galaxy. This may be 
the case for TN J0924-2201, which is undetected in SCUBA 
photometry (Reuland et al. 2004), but which according to 
our criteria has a 1.1 mm counterpart in our AzTEC map 
(§3.2) - this suggests that the AzTEC source is not the 
actual counterpart, but is instead a MMG several tens of 
kpc from the radio galaxy. Also of relevance in this context 
are results from a FWHM=2"spatial resolution SMA study 
of the z=3.79 radio galaxy 4C +60.07 (Ivison et al. 2008): 
the luminous sub-mm emission, previously assumed to be 
centred on the active galaxy's host, was resolved into two 
discrete sub-mm sources with offsets of 1"(~10 kpc) and 
4"(~30 kpc) from the position of the radio core. These two 
results provide evidence, although anecdotal in nature, that 
the millimetre/sub-mm emission associated with high-z ra- 
dio galaxies may not always be spatially coincident with the 
host galaxy. Unfortunately, with currently available facili- 
ties high resolution observations of the kind presented by 
Ivison et al. (2008) are costly and are sensitive only the very 
brightest of MMGs. The next generation of millimetre/sub- 
mm facilities, such as the 50 m Large Millimetre Telescope 
or the Atacama Large Millimetre Array, will be needed to 
resolve this issue. 

Nevertheless, we can state that all 17 of the active galax- 



ies in our sample are associated with one or more MMG, and 
thus we conclude that high-z (radio-loud) active galaxies are 
beacons for finding MMGs. 



7 SUMMARY 

We have presented 1.1 mm imaging observations of the 
7 arcmin 2 fields around 16 powerful radio galaxies at 
0.5<z<5.2 and a radio quiet quasar at z=6.3. With these 
new observations, we have more than doubled the num- 
ber of high-z radio galaxies which have been imaged at 
millimetre/sub-millimetre wavelengths. For all 17 active 
galaxies, we detect at least one associated millimetric source, 
which shows that high-z (radio-loud) active galaxies are bea- 
cons for finding millimetre/sub-millimetre galaxies at high- 
z. We identify likely millimetric counterparts for 11 radio 
galaxies, and these have 1.1 mm flux densities ranging be- 
tween <1.4 mjy (3a) and 8.6±0.9 mjy, with implied star 
formation rates of <200 to ~1300 Mq yr _ . After stacking 
our images at the position of the radio galaxies, including 
those for which we did not detect the millimetric counter- 
part, we derive an average flux density of 3.4±0.2 mjy. 

For the high-z radio galaxies, we have identified an an- 
ticorrelation between 1.1 mm flux density (or star formation 
rate) and the projected linear size of the radio source. From 
this result, we have concluded that smaller radio sources are 
associated with more powerful (obscured) starbursts. 

We have also presented images and photometry from 
new and archival Spitzer IRAC and MIPS observations. We 
find that the mid-infrared spectral energy distributions of 
the radio galaxies tend to be bluer than those of Spitzer- 
selected high-z active galaxies (Lacy et al. 2004; Martmez- 
Sansigre et al. 2008). In particular, the radio galaxies which 
we have detected at 1.1 mm have mid- infrared colours con- 
sistent with obscured star formation. We have also identified 
three ^=100 kpc arcs of 24/im sources of which three radio 
galaxies appear to be part; for the two arcs at z>2, we note 
a spatial correlation between the 24 /im and the 1.1 mm 
emission. 



ACKNOWLEDGMENTS 

We thank the referee, Carlos De Brueck, for suggestions 
which improved this paper, for making published MAMBO 
maps of 4C+41.17 and TN J1338-1942 available to us, and 
for providing us with an unpublished radio image of TN 
J2007-1316 which we used to determine its radio source posi- 
tion angle (28°) and the position of its radio core (see Table 
1). We also thank Tom Downes for his important contri- 
bution to the AzTEC data pipeline software. AH acknowl- 
edges useful discussions with Alfredo Montana and Emmaly 
Aguilar, and financial support from both CONACyT and 
INAOE. This research was partially supported by CONA- 
CyT project 39953F. This work has also been supported in- 
part by NSF Grants 0907952 and 0838222. The James Clerk 
Maxwell Telescope is operated by The Joint Astronomy Cen- 
tre on behalf of the Science and Technology Facilities Coun- 
cil of the United Kingdom, the Netherlands Organisation 
for Scientific Research, and the National Research Council 
of Canada. The ASTE project is driven by the Nobeyama 



16 A. Humphrey et al. 



Radio Observatory (NRO), a branch of the National As- 
tronomical Observatory of Japan (NAOJ), in collaboration 
with the University of Chile and Japanese institutions in- 
cluding the University of Tokyo, Nagoya University, Osaka 
Prefecture University, Ibaraki University, and Hokkaido Uni- 
versity. 



REFERENCES 

Archibald E. N., Dunlop J. S., Hughes D. H., Rawlings S., 
Eales S. A., Ivison R. J., 2001, MNRAS, 323, 417 
Athreya R. M., Kapahi V. K., McCarthy P. J., van Breugel 
W., 1998, A&A, 329, 809 (At98) 
Austermann J. E., et al., 2010, MNRAS, 401, 160 
Becker R. H., White R. L., Edwards A. L., 1991, ApJS, 75, 
1 (Be91) 

Benford D. J., 1999, Thesis, California Institute of Technol- 
ogy (Be99) 

Bertin E. & Arnouts S., 1996, A&AS, 117, 393 
Best P. N., Bailer D. M., Longair M. S., Riley J. M., 1995, 
MNRAS, 275, 1171 

Best P. N., Longair M. S., Rottgering H. J. A., 1996, MN- 
RAS, 280, L9 

Best P. N., Rottgering H. J. A., Lehnert M. D., 1999, MN- 
RAS, 310, 223 (Be99) 

Best P. N., Rottgering H. J. A., Longair M. S., 2000, MN- 
RAS, 311, 23 

Binette L., Wang J., Villar-Martin M., Martin P. C, Magris 
C. C, 1993, ApJ, 414, 535 

Black A. R. S., Baum S. A., Leahy J. P., Perley R. A., Riley 
J. M., Scheuer P. A. C, 1992, MNRAS, 256, 186 
Blundell K. M., 1996, MNRAS, 283, 538 
Bolton J. C, Shimmins A. J., Wall J. V., 1975, AuJPA, 34, 
1 (Bo75) 

Bornancini C. C, De Breuck O, de Vries W., Croft S., van 
Breugel W., Rottgering H., Minniti D., 2007, MNRAS, 378, 
551 (Bo07) 

Broderick J. W., De Breuck O, Hunstead R. W., Seymour 
N., 2007, MNRAS, 375, 1059 (Br07) 

Burgess A. M., Hunstead R. W., 2006, AJ, 131, 114 (Bu06) 
Carilli C. L., Roettgering H. J. A., van Ojik R., Miley G. K., 
van Breugel W. J. M., 1997, ApJS, 109, 1 (Ca97) 
Carilli C. L., Wang R., van Hoven M. B., Dwarakanath K., 
Chengalur J. N., Wyithe S., 2007, AJ, 133, 2841 (Ca07) 
Cimatti A., di Serego Alighieri S., Vernet J., Cohen M., Fos- 
bury R. A. E., 1998, ApJ, 499, L21 (Ci98) 
Chambers K. O, Miley G. K., van Breugel W., 1987, Natur, 
329, 604 

Chambers K. O, Miley G. K., van Breugel W. J. M., 1990, 
ApJ, 363, 21 (Ch90) 

Chambers K. C, Miley G. K., van Breugel W. J. M., Bre- 
mer M. A. R., Huang J.-S., Trentham N. A., 1996, ApJS, 
106, 247 (Ch96) 

Chen, Y., Lowenthal, J. D., Yun, M. S., 2010, ApJ, 712, 1385 
Chini R., Kriigel E., Kreysa E., 1986, A&A, 167, 315 
Cohen A. S., Lane W. M., Cotton W. D., Kassim N. E., 
Lazio T. J. W., Perley R. A., Condon J. J., Erickson W. C, 
2007, AJ, 134, 1245 (Co07) 

Condon J. J., Cotton W. D., Greisen E. W., Yin Q. F., Per- 
ley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 1693 
(Co98) 



De Breuck C, van Breugel W., Minniti D., Miley G., 

Rottgering H., Stanford S. A., Carilli C, 1999, A&A, 352, 

L51 (DeBr99) 

De Breuck C, Rottgering H., Miley G., van Breugel W., 

Best P., 2000a, A&A, 362, 519 

De Breuck C, van Breugel W., Rottgering H. J. A., Miley 

G., 2000b, A&AS, 143, 303 (DeBrOO) 

De Breuck C, et al., 2001, AJ, 121, 1241 (DeBrOl) 

De Breuck C, van Breugel W., Stanford S. A., Rottgering 

H., Miley G., Stern D., 2002, AJ, 123, 637 (DeBr02) 

De Breuck C, et al., 2004, A&A, 424, 1 (DeBr04) 

De Breuck C, et al., 2010, ApJ, 725, 36 

Douglas J. N., Bash F. N., Bozyan F. A., Torrence G. W., 

Wolfe C, 1996, AJ, 111, 1945 (Do96) 

Downes A. J. B., Peacock J. A., Savage A., Carrie D. R., 

1986, MNRAS, 218, 31 

Downes D., Solomon P. M., Sanders D. B., Evans A. S., 

1996, A&A, 313, 91 

Dunlop J. S., Hughes D. H., Rawlings S., Eales S. A., Ward 

M. J., 1994, Natur, 370, 347 

Dunlop, J., Peacock, J., Spinrad, H., Dey, A., Jimenez, R., 

Stern, D., & Windhorst, R., 1996, nat, 381, 581 

Efstathiou A., Rowan- Robinson M., Siebenmorgan R., 2000, 

MNRAS, 313, 734 

Fan X., et al., 2001, AJ, 122, 2833 

Fanti C, Fanti R., Ficarra A., Padrielli L., 1974, A&AS, 18, 

147 (Fa74) 

Fazio G. G. et al. 2004, ApJSS, 154, 10 

Ficarra A., Grueff G., Tomassetti G., 1985, A&AS, 59, 255 

(Fi85) 

Foley R. J., et al., 2011, arXiv.1101.1286 

Gopal-Krishna, Ledoux C, Melnick J., Giraud E., Kulkarni 

V., Altieri B., 2005, A&A, 436, 457 (Go05) 

Gower J. F. R., Scott P. F., Wills D., 1967, MmRAS, 71, 49 

(Go67) 

Griffith M., Langston G., Heflin M., Conner S., Lehar J., 

Burke B., 1990, ApJS, 74, 129 (Gr90) 

Gregory P. C, Condon J. J., 1991, ApJS, 75, 1011 (Gr91) 

Gregory P. C, Vavasour J. D., Scott W. K., Condon J. J., 

1994, ApJS, 90, 173 (Gr94) 

Greve T. R., Stern D., Ivison R. J., De Breuck C, Kovacs 

A., Bcrtoldi F., 2007, MNRAS, 382, 48 (Gr07) 

Griffith M. R., Wright A. E., Burke B. F., Ekers R. D., 1994, 

ApJS, 90, 179 (Gr94) 

Griffith M. R., Wright A. E., Burke B. F., Ekers R. D., 1995, 

ApJS, 97, 347 (Gr95) 

Gutermuth R. A., Myers P. C, Megeath S. T., Allen L. E., 

Pipher J. L., Muzerolle J., Porras A., Winston E., Fazio G., 

2008, ApJ, 674, 336 

Hales S. E. G., Baldwin J. E., Warner P. J., 1993, MNRAS, 

263, 25 (Ha93) 

Hatch N. A., Overzier R. A., Kurk J. D., Miley G. K., 

Rottgering H. J. A., Zirm A. W., 2009, MNRAS, 395, 114 

Healey S. E., Romani R. W., Taylor G. B., Sadler E. M., 

Ricci R., Murphy T., Ulvestad J. S., Winn J. N., 2007, ApJS, 

171, 61 (He07) 

High F. W. et al., 2010, ApJ, 723, 1736 

Hildebrand R. H., 1983, QJRAS, 24, 267 

Humphrey A., Villar-Martin M., Fosbury R., Vernet J., di 

Serego Alighieri S., 2006, MNRAS, 369, 1103 

Humphrey A., Villar-Martin M., Fosbury R., Binette L., 

Vernet J., De Breuck C, di Serego Alighieri S., 2007, MN- 



17 



RAS, 375, 705 

Humphrey A., Villar-Martin M., Vernet J., Fosbury R., di 
Serego Alighieri S., Binette L., 2008, MNRAS, 383, 11 
Inskip K. J., Best P. N., Rottgering H. J. A., Rawlings S., 
Cotter G., Longair M. S., 2002, MNRAS, 337, 1407 
Isobe T., Feigelson E. D., Nelson P. I., 1986, ApJ, 306, 490 
Ivison R. J., Dunlop J. S., Smail Ian, Dey Arjun, Liu Michael 
C, Graham J. R., 2000, ApJ, 542, 27 
Ivison R. J., et al., 2007, MNRAS, 380, 199 
Ivison R. J., et al., 2008, MNRAS, 390, 1117 
Jeyakumar S., Saikia D. J., Pramesh Rao A., Balasubrama- 
nian V., 2000, A&A, 362, 27 (JeOO) 
Kim S., et al., 2009, ApJ, 695, 809 

Knopp G. P., Chambers K. C, 1997, ApJS, 109, 367 (Kn97) 
Kuehr H., Witzel A., Pauliny-Toth I. I. K., Nauber U., 1981, 
A&AS, 45, 367 (Ku81) 
Lacy M., et al., 2004, ApJS, 154, 166 

Large M. I., Mills B. Y., Little A. G., Crawford D. F., Sut- 
ton J. M., 1981, MNRAS, 194, 693 (La81) 
Liu X., Cui L., Luo W.-F., Shi W.-Z., Song H.-G., 2007, 
A&A, 470, 97 (Li07) 

Martfnez-Sansigre A., Lacy M., Sajina A., Rawlings S., 2008, 
ApJ, 674, 676 

Mauch T., Murphy T., Buttery H. J., Curran J., Hunstead 
R. W., Piestrzynski B., Robertson J. G., Sadler E. M., 2003, 
MNRAS, 342, 1117 (Ma03) 

McCarthy P. J., van Breugel W., Spinrad H., Djorgovski S., 
1987, ApJ, 321, L29 

McCarthy P. J., Spinrad H., van Breugel W., 1995, ApJS, 
99, 27 

McCarthy P. J., Kapahi V. K., van Breugel W., Persson 
S. E., AthreyaR., Subrahmanya C. R., 1996, ApJS, 107, 19 
(Mc96) 

McLure, R. J., Kukula, M. J., Dunlop, J. S., Baum, S. A., 
O'Dea, C. P., & Hughes, D. H., 1999, MNRAS, 308, 377 
Menanteau F. et al., 2010, ApJ, 723, 1523 
Miley G. K., Chambers K. C, van Breugel W. J. M., Mac- 
chetto F., 1992, ApJ, 401, L69 (Mi92) 
Miley G., De Breuck C, 2008, A&ARv, 15, 67 
Murgia M., Fanti C, Fanti R., Gregorini L., Klein U., Mack 
K.-H., Vigotti M., 1999, A&A, 345, 769 
Muxlow T. W. B., Pelletier G., Roland J., 1988, A&A, 206, 
237 

Nesvadba N. P. H., et al., 2006, ApJ, 650, 661 
Nesvadba N. P. H., et al., 2009, MNRAS, 395, L16 
Overzier R. A., Rottgering H. J. A., Kurk J. D., De Breuck 
C, 2001, A&A, 367, L5 (OvOl) 

Pentericci L., Roettgering H. J. A., Miley G. K., Carilli C. L., 
McCarthy P., 1997, A&A, 326, 580 (Pe97) 
Pentericci L., Van Reeven W., Carilli C. L., Rottgering 
H. J. A., Miley G. K., 2000, A&AS, 145, 121 (PeOO) 
Pentericci L., McCarthy P. J., Rottgering H. J. A., Miley 
G. K., van Breugel W. J. M., Fosbury R., 2001, ApJS, 135, 
63 

Pilkington J. D. H., Scott P. F., 1965, MmRAS, 69, 183 
(Pi65) 

Priddey R. S., Isaak K. G., McMahon R. C, Robson E. I., 
Pearson C. P., 2003, MNRAS, 344, L74 
Priddey R. S., Ivison R. J., Isaak K. C, 2008, MNRAS, 383, 
289 

Reuland M., et al., 2003a, ApJ, 592, 755 (Re03) 
Reuland Michiel, van Breugel Wil, Rottgering Huub, de 



Vries Wim, De Breuck Carlos, Stern Daniel, 2003b, ApJ, 

582L, 71 

Reuland M., Rottgering H., van Breugel W., De Breuck C, 

2004, MNRAS, 353, 377 (Re04) 

Reuland Michiel, 2005, PhD Thesis, Leiden Observatory, 

Leiden University 

Rickett B. J., Lazio T. J. W., Ghigo F. D., 2006, ApJS, 165, 

439 (Ri06a) 

Ricci R., Prandoni I., Gruppioni C, Sault R. J., de Zotti 

C, 2006, A&A, 445, 465 (Ri06b) 

Rieke G. H. et al., 2004, ApJSS, 154, 25 

Rottgering, H. J. A., van Ojik, R., Miley, G. K., Chambers, 

K. C, van Breugel, W. J. M., & de Koff, S., 1997, A&A, 

326, 505 (R697) 

Scott K. S., et al., 2010, MNRAS, 405, 2260 

Seymour N., et al., 2007, ApJS, 171, 353 

Slee O. B., 1995, AuJPh, 48, 143 (S195) 

Solorzano-Inarrea C, Best P. N., Rottgering H. J. A., 

Cimatti A., 2004, MNRAS, 351, 997 

Stanghellini C, O'Dea C. P., Dallacasa D., Baum S. A., 

Fanti R., Fanti C, 1998, A&AS, 131, 303 (St98) 

Stanghellini C, O'Dea C. P., Dallacasa D., Cassaro P., 

Baum S. A., Fanti R., Fanti C, 2005, A&A, 443, 891 (St05) 

Stern D., et al., 2005, ApJ, 631, 163 

Stevens J. A., et al., 2003, Natur, 425, 264 

Stevens J. A., Jarvis M. J., Coppin K. E. K., Page M. J., 

Greve T. R., Carrera F. J., Ivison R. J., 2010, MNRAS, 405, 

2623 

Stiavelli M., et al., 2005, ApJ, 622, LI 

van Breugel W., Miley G., Heckman T., Butcher H., Bridle 

A., 1985, ApJ, 290, 496 

van Breugel W., De Breuck C, Stanford S. A., Stern D., 

Rottgering H., Miley C, 1999, ApJ, 518, L61 (vB99) 

van Ojik R., Roettgering H. J. A., Miley G. K., Hunstead 

R. W., 1997, A&A, 317, 358 

Venemans B. P., et al., 2007, A&A, 461, 823 (Ve07) 

Vernet J., Fosbury R. A. E., Villar-Martin M., Cohen M. H., 

Cimatti A., di Serego Alighieri S., Goodrich R. W., 2001, 

A&A, 366, 7 

Vieira J. D., et al., 2010, ApJ, 719, 763 

Villar-Martin M., Fosbury R. A. E., Binette L., Tadhunter 

C. N., Rocca-Volmerange B., 1999, A&A, 351, 47 (Vi99) 

Villar-Martin M., Vernet J., di Serego Alighieri S., Fosbury 

R., Humphrey A., Pentericci L., 2003, MNRAS, 346, 273 

Villar-Martin M., et al., 2006, MNRAS, 366, LI 

Villar-Martin M., Humphrey A., De Breuck C, Fosbury R., 

Binette L., Vernet J., 2007, MNRAS, 375, 1299 

Viner M. R., Erickson W. C, 1975, AJ, 80, 931 (Vi75) 

Werner M. W. et al., 2004, ApJSS, 154, 1 

White R. L., Becker R. H., 1992, ApJS, 79, 331 (Wh92) 

Willott C. J., Rawlings S., Archibald E. N., Dunlop J. S., 

2002, MNRAS, 331, 435 

Willott C. J., Delfosse X., Forveille T., Delorme P., Gwyn 

S. D. J., 2005, ApJ, 633, 630 

Wilson G. W., et al., 2008, MNRAS, 386, 807 

Wright A., Otrupcek R., 1990, Parkes Catalogue, Australia 

Telescope National Facility (Wr90) 

Wright A. E., Griffith M. R., Burke B. F., Ekers R. D., 1994, 

ApJS, 91, 111 (Wr94) 

Xiang L., Dallacasa D., Cassaro P., Jiang D., Reynolds C, 

2005, A&A, 434, 123 (Xi05) 

Xiang L., Reynolds C, Strom R. C, Dallacasa D., 2006, 



A. Humphrey et al. 



A&A, 454, 729 (Xi06) 

Yun M. S., et al., 2008, MNRAS, 389,333 

Yun M. S., et al., 2011, MNRAS, submitted 



